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Exact  Theory  of  One-Phonon  Amplitudes  in  Diffractive  Atom 
Surface  Scattering 

Michael  E.  Flatte 


This  paper  is  an  extension  of  a  previous  paper,  Exact  Theory  of 
One-Phonon  Amplitudes  in  Atom-Surface  Scattering,  by  Michael 
E.  Flatte  and  Walter  Kohn.  That  paper  derives  an  expression  for  the 
probabihty  of  creating  or  annihilating  one  long-wavelength  phonon  in 
an  atom-surface  scattering  event  for  a  system  with  negligible  diffraction. 
The  expression  is  independent  of  the  atom-surface  potential  and  of  force 
constant  changes  near  the  surface.  It  is  exact  if  the  inelastic  scattering  is 
weak.  This  paper  derives  a  result  with  the  same  qualitative  features  for 
systems  with  arbitrary  diffraction. 

Determining  the  atom-surface  potential  is  a  complicated  problem 
which  is  very  important  for  such  fields  as  corrosion,  catalysis  and  surface 
diffusion.  Atom-surface  scattering  experiments  are  good  probes  of  this 
potential.  It  is  useful  for  an  experimenter  to  know  that  some  phonons 
will  be  ineffective  in  probing  this  potential.  Our  paper  proves  that  this 
category  includes  all  long-wavelength  phonons  when  the  total  inelastic 
scattering  is  weak. 

Furthermore,  an  exact  expression  for  the  scattering  probability  which 
is  independent  of  the  atom-surface  potential  may  be  useful  for  calibrating 
other  experiments. 
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Exact  Theory  of  One-Phonon  Amplitudes  in  Diffractive  Atom- 
Surface  Scattering 

Michael  E.  Flatte 

Department  of  Physics,  University  of  California,  Santa  Barbara,  CA 
93106 

(To  be  published  in  the  proceedings  of  the  Third  International  Con¬ 
ference  on  the  Structure  of  Surfaces  (ICSOS-3)  ) 

Abstract  I  derive  an  expression  for  the  probabihty  of  creating  or 
annihilating  one  long-wavelength  surface  or  bulk  phonon  during  a 
scattering  event  in  a  system  with  diffraction  scattering.  This  ex¬ 
pression  depends  on  the  hulk  elastic  constants  but  is  independent 
of  force  constant  changes  near  the  surface  and  of  the  details  of  the 
atom-target  potential.  If  the  inelastic  scattering  is  weak  (as  defined 
in  the  text),  the  expression  is  exact. 

In  a  system  where  the  only  elastic  scattering  is  specular,  the  bulk 
elastic  constants  of  the  target  and  the  quantum  numbers  of  the  decou¬ 
pled  atom-target  system  entirely  deternaine  the  probability  for  the  weak 
inelastic  creation  (or  annihilation)  of  a  long-wavelength  phonon  [  1].  A 
heuristic  explanation  follows  of  why  an  expression  with  similar  qualities 
should  also  exist  for  systems  with  non-specular  elastic  scattering.  In  par¬ 
ticular,  the  expression  is  independent  of  (1)  deviations  of  force  constants 
from  bulk  values  near  the  surface  of  a  target  uncoupled  to  atoms  and  (2) 
the  fimctional  form  of  the  atom-target  interaction  potential  F(r;  {{?,}). 
where  r  is  the  position  of  the  colliding  atom  and  Ui  the  displacement  of 
the  z’th  target  atom. 

Displacements  in  the  surface  region  associated  with  a  phonon  of  small 
w'ave-vector  q  are  quasi-rigid  for  {qo,)~^  layers  perpendicular  to  the  sur¬ 
face,  where  a  is  the  lattice  constant.  If  ga  -C  1,  force  constant  changes  in 
a  few  layers  near  the  surface  are  irrelevant  to  the  displacements  associated 
with  this  phonon. 

The  functional  form  of  V  is  irrelev2Lnt  because  to  lowest  order  in  qa. 
the  first-order  expression  for  the  probability  is  proportional  to 


Here  V{f^  =  U(r;  {0})  and  the  |’^'o)’s  are  wave  functions  of  the  atom-rigid 
tcirget  system  with  the  same  crystal  momentum  parallel  to  the  surface 


(K)  and  energy  (E).  When  there  is  only  elastic  scattering,  the  case  with 
a  rigid  target,  for  a  given  K  and  E  there  is  a  discrete  set  of  asymptotic 
perpendicular  momenta  the  atom  may  have.  These  channels,  incoming 
if  <  0,  outgoing  if  >  0,  are  labelled  by  surface  reciprocal  lattice 
vectors  G.  The  atom  involved  in  an  inelastic  scattering  process  where 
ga  <C  1  Will  emerge  near  one  of  these  channels.  The  probability  that  the 
atom  emerges  near  channel  Go  is  calculated  by  choosing  the  amplitude 
of  all  incoming  channels  except  0  to  be  0  for  and  the  amplitude 

of  all  outgoing  channels  except  Go  to  be  0  for  1^©^^).  When  both  wave 
functions  are  normalized  to  unit  incoming  flux  in  channel  0,  the  integral 
represents  the  total  force  applied  in  transferring  the  atom  from  incoming 
channel  0  to  outgoing  channel  Go-  This  is  rigorously  given  by  the  change 
in  perpendicular  momentum  per  unit  time  and  area,  i.e.,  +  k^). 

The  results  axe  exact  in  the  limit  where  the  following  conditions  are 
satisfied: 

kBT<hujD  (2d)  ga  <  1  (2e), 

where  c  is  a  characteristic  velocity  of  a  phonon  in  the  target,  h  is  the 
range  of  the  atom-target  potential  and  Vo  is  its  depth,  m  is  the  mass 
of  the  colliding  atom  and  Ea  is  its  incident  energy.  M  is  the  mass  of  a 
target  atom,  ljd  is  the  Debye  frequency,  and  T  is  the  temperature  of  the 
target. 

The  condition  (2a)  allows  one  to  stop  the  expansion  of  V  after  the 
first  order  term 

\ 

F(r-;  {5.})  =  V{f)  +  Si  ■  {0})  (3) 

Condition  (2b)  ensures  that  in  calculating  the  inelastic  scattering,  the 
last  term  in  (3)  can  be  treated  to  lowest  order. 

Derivation  of  the  Probability  for  Production  of  Long- Wavelength 
Phonons 

This  derivation  requires  the  use  of  a  variational  principle  developed  for 
multichannel  scattering  in  nuclear  collisions  [  2].  The  following  assump¬ 
tions  simplify  this  presentation:  (1)  the  atom  only  interacts  with  the 
uppermost  layer  of  target  atoms  and  (2)  the  units  are  chosen  so  that  the 
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target’s  volume  and  the  area  of  its  interacting  surface  axe  1 .  The  results 
remain  unchanged  if  these  assumptions  are  dropped. 

To  first  order  in  {wt},  the  Hamiltonian  consists  of  an  unperturbed 
Hamiltonian  Ho->  which  ignores  phonon  coupling  to  the  colliding  atom, 
and  Hi,  which  couples  the  atom  linearly  to  the  phonons. 


H  =  ifportioi.  +  Hiarget  +  V(f,  {u.})  =  H„  +  Hi  +  0(«2)  (4) 


where 


^  ^  +  V(r-)  (5a)  (56) 


i.a 


Here  a  denotes  a  cartesian  coordinate  and  is  the  occupa¬ 

tion  of  phonon  mode  q,  v,  where  v  is  the  phonon’s  branch. 

Let  and  be  two  wave  functions  for  the  combined  target- 

atom  system  with  the  same  E  and  K.  Their  asymptotic  form  is 


{n},(5 


^  G,{nr 


(6) 


where  the  sum  is  carried  out  over  all  surface  reciprocal  lattice  vectors  G 
and  phonon  occupations  {n},  R  is  the  atom’s  displacement  parallel  to 
the  surface,  and 


2m 


/?,„}  =  K-Y,  +  Gf  =  -^{E  -  Y, 


,  and  Then  the  equation 

G,{n}  G,{n}  ^ 


(7) 

Choose  two  trial  functions  and  such  that  = 

>(1)  _  d(1) 


/(«<*>|r)(if -EXrl^f  (8) 

vW  Art 


<5,{n} 


where  6B^ \  ,  =  B^3\  ,  -  , ,  is  correct  to  first  order  in  the  errors 

G,{n}  G,{n},t  G,{n}’ 

of  the  wave  functions. 
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This  variational  principle  is  of  value  because  the  integral  can  be 
evaluated  analytically  to  lowest  order  in  qa  for  trial  functions  which  are 
eigenfunctions  of  Ho>  These  wavefunctions  have  the  form  [  3] 


(9) 


The  set  of  all  <f>(^{z)  for  a  given  E  and  {n}  satisfy  the  following  coupled 

_  ^ 

equations,  where  Ea  —  E  -  and  =  V{f): 


2m  “  2m  ^ 


(10) 

Let  o)  ^  eigenfunction  of  Hq  with  energy  E,  crystal  mo¬ 
mentum  K  and  l{n2})  =  1  •  •  •  . . .)  for  which  and  all 

other  j4’s  are  0.  Let  o)  be  an  eigenfunction  of  Ho  with  eigenvalues 

E,  K  and  |{ni})  =  ] . . .  -1-  1 . . .)  for  which  and 

is  the  only  nonzero  B. 

For  these  trial  functions,  the  integral  is  }  ol^^i  }  o)^^- 

Expanding  the  Ui  in  normal  modes  and  rewriting  yields 


q,r/,Q 


dV 

f  dUi,a 


{r}e 


-f  c.c. 


(11) 


The  second  exponential  can  be  ignored  if  the  wavelength  of  the  phonon 
is  much  larger  than  the  range  of  the  coupling  potential,  similaj  to  a 
and  6,  since  unless  q{Ri  —  R)  1  (dV/dui)  will  vanish.  Furthermore, 
'^idV/dui^a{r)  =  —dV{f)fdra  since  V(r;  {^i})  depends  only  on  f  —  Ui. 

t 

Finally,  defining  +  ^-qA  yields 


(12) 


q,u,a 


Transforming  the  resulting  integral  with  (5a)  and  Green’s  Theorem 
yields 

d  . 


(nj,,  +  ^  j^dS^ 


G,G',a 


where  h  is  the  unit  normal  to  the  surface  S  which  bounds  V.  Because 
of  the  smixl^  difference  in  atomic  energy  and  parallel  crystal  momentum 
between  and  may  be  replaced  by 

The  vanishing  of  the  resulting  integrals  for  Uxid/dx)  and  ny{d/dyj  and 
evaluation  of  the  remainder  yields 


+  1)^ 


2m 


E 


<5  dr  adz  dz  dr  a 


(14) 


^  ^=oo 


Since  (14)  is 

zero  if  a  denotes  x  or  and  when  a  denotes  2  it  is 

^Go)- 

The  differential  probability  of  scattering  an  atom  into  final  angle  fl 
with  energy  Ea-,  per  unit  incident  perpendicular  flux  is 


(J?P 


dSldEa  (15) 


where  k  is  the  atom’s  final  momentum,  p{Ea)  is  the  density  per  unit 
energy  of  atomic  states  with  momentum  in  final  angle  Q,  and  g(hu;^y,  Q) 
is  the  density  per  unit  energy  of  phonons  with  frequency  uj  =  and 
surface  wavevector  Q. 

For  a  target  which  is  isotropic  in  the  long  wavelength  limit, 
and  g{fiu;^y,Q)  are  well  known,  depend  only  on  bulk  elastic  constants, 
density,  u)  and  Q,  and  are  given  in  sample  form  in  Ref.  1.  For  non¬ 
isotropic  targets,  U^^u,z  and  g{huj^y,Q)  depend  in  a  more  complicated 
way  onJ;he  bulk  elastic  constants,  density.  Miller  indices  of  the  surface, 
u!  and  Q. 
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